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I. INTRODUCTION 

BCS theory predicts that the phase transition from the 
normal state to the superconducting state in s-wave su- 
perconductors is continuous or second order. However, 
in 1974 Halperin, Lubensky, and Ma 1 argued, based on 
a 4 — e expansion, that the coupling between the super- 
conducting order parameter and the electromagnetic vec- 
tor potential drives the transition first order. This con- 
clusion is inevitable for extreme type-I superconductors 
where fluctuations of the order parameter are negligible 
and the vector potential can be integrated out exactly. 
The mechanism in this case is known as a fluctuation- 
induced first order transition, and it is analogous to the 
spontaneous mass generation known in particle physics as 
the Coleman- Weinberg mechanism* 2 . When order param- 
eter fluctuations cannot be neglected, and especially for 
type-II superconductors, the problem cannot be solved 
exactly. The authors of Ref. [l| generalized the problem 
by considering an n/2-dimensional complex order param- 
eter and conducting a renormalization-group (RG) anal- 
ysis in d — 4 — e dimensions. The physical case of interest 
is n — 2 and e = 1. To first order in e they found that 
a RG fixed point corresponding to a continuous phase 
transition exists only for n > 365.9, which suggests that 
for physical parameter values the transition is first order 
even in the type-II case. They corroborated this conclu- 
sion by performing a large-n expansion for fixed d = 3. 
To first order in 1/n, the critical exponent v is positive 
only for n > 9.72, which again strongly suggests that the 
transition in the physical case n = 2 is first order. For 
superconductors, due to the large correlation length and 
the correspondingly small size of the fluctuations, the size 
of the first order discontinuity is too small to be observ- 
able. For the analogous smectic-A to nematic transition 
in liquid crystals, on the other hand, it was predicted to 
be much larger, and, indeed, experimentally observable. 
In contrast to this theoretical prediction, experiments on 
liquid crystals showed, and continue to show to this day, 
a clear second order transition^ This prompted a re- 
examination of the theory by Dasgupta and Halperin. 4 
Using Monte Carlo data and duality arguments, these 
authors argued that a strongly type-II superconductor in 
d = 3 should show a second order transition after all. A 



Monte-Carlo study of the intermediate region suggested 
that the transition is first order in the strongly type-I re- 
gion and continuous in the strongly type-II region, with 
a tricritical point separating the two regimes^ Why the 
continuous transition does not show as a critical fixed 
point in a 4 — e expansion is not quite clear. One pos- 
sible explanation is that the fixed point is not perturba- 
tively accessible. Another is that the critical value of n 
in the generalized model, which is close to 366 near d = 4 
where the e-expansion is controlled, decreases to a value 
less than 2 as the dimension is lowered to the physical 
value d = 3. Herbut and Tesanovio^ have shown that 
the critical value of n, above which there is a second order 
transition, decreases rapidly with increasing e, and that 
a RG analysis in a fixed dimension d — 3 to one-loop or- 
der does yield a critical fixed point for systems that are 
sufficiently strongly type II. 

Recently there has been substantial interest in 
unconventional superconductivity. In particular, 
Sr2RuC>4 has emerged as a convincing case of p-wave 
superconductivity^ and UGe2 is another candidate. 10 
This raises the question whether for such systems the 
fluctuation-induced first order mechanism is also appli- 
cable, or whether the additional order parameter degrees 
of freedom allow for a critical RG fixed point, signaliz- 
ing a second order transition, even though no such fixed 
point is found in the s-wave case. Here we investigate 
this problem. By conducting an analysis for p-wave su- 
perconductors analogous to the one of Ref. [l] we find that 
there is no critical fixed point, as in the s-wave case. This 
analysis thus also suggests a first order transition, as in 
the s-wave case, although the restrictions are somewhat 
less stringent. Presumably, the same reservations regard- 
ing non-perturbative fixed points that are suspected to 
be relevant for the s-wave case apply here as well. 

This paper is organized as follows. In Sec. HI] we de- 
fine our model and derive the mean-field phase diagram. 
In Sec. Mil we determine the nature of the phase transi- 
tion. We do so first in a renormalized mean-field approxi- 
mation that neglects fluctuations of the superconducting 
order parameter. We then take such fluctuations into 
account, first by means of a RG analysis in d = 4 — e 
dimensions, and then by means of a 1/n-expansion. In 
Sec. HVI we summarize our results. 
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II. MODEL 



III. NATURE OF THE PHASE TRANSITION 



Let us consider a Landau-Ginzburg- Wilson (LGW) 
functional appropriate for describing spin-triplet super- 
conducting order. The superconducting order parame- 
ter is conveniently written as a matrix in spin spaced 
A a/ 3(fc) = X^Ui d »( k ) ( f7 M* cr 2) Q/3 - Here cri j2 ,3 are the 
Pauli matrices, k is a wave vector, and the are the 
components of a complex 3- vector d(k). p-wave sym- 
metry implies d fl (k) — Y^j=i d nj^j-> with k a unit wave 
vector. The tensor field d fl j(x) is the general order pa- 
rameter for a spin-triplet p-wave superconductor and it 
allows for a very rich phenomenology. For dcfiniteness we 
will constrain our discussion to a simplified order param- 
eter describing the so-called /3-state^ which has been 
proposed to be an appropriate description of UGe2 It 
is given by a tensor product d — tp ® <fi of a complex 
vector i/> in spin space and a real unit vector cj> in or- 
bital space. The ground state is given by tp — A (l, i, 0), 
4> = (0, 0, 1). In a weak-coupling approximation that ne- 
glects terms of higher than bilinear order in ip 2 , <p 2 , and 
V 2 the action depends only on ip, 



dx 



* I 2 



t\ip\ 2 + c\Dxp\ 2 + u\ip\ 4 + v\xp xf| 



(2.1) 



Here A is the vector potential, D = V iqA is the 
gauge invariant gradient with q = 2e the Cooper pair 
charge (we use units such that Planck's constant and the 
speed of light are unity), and \Dip\ 2 = (Diip a )(D*ip^) 
with summations over i and a implied, /i is the normal- 
state magnetic permeability, and t, c, u, and v are the 
parameters of the LGW functional. The fields ip and A 
are understood to be functions of the position x. 

For later reference we now generalize the vector ip from 
a complex 3- vector to a complex m- vector with compo- 
nents tp a , so that the total number of order-parameter 
degrees of freedom is n = 2m. In order to generalize the 
term with coupling constant v we use of the following 
identity for 3- vectors, 



\lp X VT = 0aCW? - 4>a1pa^0lp} 



(2.2) 



and notice that the right-hand side is well defined for a 
complex m-vector. Our generalized action now reads 



S = 



dx 



t ip a ?p* a + c (AVOWO - « ^o^li>* 



+ (u + v) ipaiplippipp + ^— e l j k {djA k )e am (diA m ) 

(2.3) 

with a, (3 — l,...m; ... — 1,2,3; and summation 
over repeated indices implied. In addition to the gener- 
alization of the order parameter to an m-vector we will 
also consider the system in a spatial dimension d close to 
d = 4. The physical case of interest is m = d = 3. 



A. Mean-field approximation 

The simplest possible approximation ignores both the 
fluctuations of the order parameter field if) and the elec- 
tromagnetic fluctuations described by the vector poten- 
tial A. The order parameter is then a constant, ip(x) = 
ip, and the free energy density / reduces to 



f = t\iP\ 2 + u\iP\' l + v\iPxiP*\ 



*|2 



(3.1) 



In order to determine the phase diagram we parame- 
terize the order parameter as follows^ 



(3.2) 



Here ipo is a real-valued amplitude, n and rh are indepen- 
dent real unit vectors, and is a phase angle. The free 
energy density can then be written 

f = tipl + u + v </>o (l - (ft • m) 2 ) sin 2 20 . (3.3) 

We now need to distinguish between two cases. 

Case 1: v > 0. The free energy is minimized by h = m, 
and ipo = —t/2u. The condition u > must be fulfilled 
for the system to be stable. 

Case 2: v < 0. The free energy is minimized by n _L rh 
and 4> = 7r/4, and ipo — —t/2(u + v). The condition 
u + v > must be satisfied for the system to be stable. 

The first case implies xpxxp* — 0. This is referred to as 
the unitary phase. In the second case, ip X ip* ^ 0, which 
is referred to as the non-unitary phase. In either case, 
mean-field theory predicts a continuous phase transition 
from the disordered phase to an ordered phase at t = 0. 
The mean-field phase diagram in the u-v plane is shown 
in Fig. [TJ 



B. Renormalized mean-field theory 

A better approximation is to still treat the order pa- 
rameter as a constant, ip(x) = but to keep the elec- 
tromagnetic fluctuations. The part of the action that 
depends on the vector potential then takes the form 



S A = ^-Jdx [k 2 A 2 (x) + (V x A(x))] 



where 



2W2 



(3.4a) 



(3.4b) 



is the inverse London penetration depth. Since A enters 
Sa only quadratically, it can be integrated out exactly^ 
and the technical development is identical to the s-wave 
caseJ*^ The result for the leading terms in powers of 
\ip\ 2 in d = 3 is 

/ = t \*P\ 2 + u H 4 - w (M 2 ) 3/2 + v \iP x ?p*\ 2 . (3.5) 
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FIG. 1: Mean-field phase diagram of a p-wave superconduc- 
tor as described by Eq. (|2.ip . See the text for additional 
information. 



Here w oc sj \xq 2 is a positive coupling constant whose 
presence drives the transition into either of the ordered 
phases first order. 

There are several interesting aspects of this result. 
First, the additional term in the mean-field free energy, 
with coupling constant w, is not analytic in \ip\ . This 
is a result of integrating out the vector potential, which 
is a soft or massless fluctuation. Second, the resulting 
first-order transition is an example of what is known as 
the Coleman- Weinberg mechanism in particle physics, 2 
or a fluctuation-induced first-order transition in statisti- 
cal mechanics X 

Let us discuss the validity of the renormalized mean- 
field theory. The length scale given by the London pen- 
etration depth A = fc^ 1 needs to be compared with 
the second length scale that characterizes the action, 
Eq. ()2.1|) . which is the superconducting coherence length 
£ = ycJWV The ratio n — A/£ is the Landau-Ginzburg 
parameter. For n — > 0, order parameter fluctuations are 
negligible (this is the limit of an extreme type-I supercon- 
ductor), and the renormalized mean-field theory become 
exact. For nonzero values of k the fluctuations of the 
order parameter cannot be neglected, and the question 
arises whether or not they change the first-order nature 
of the transition. We will investigate this question next 
by means of two different technical approaches. 



Eq. (|2.3p . For the ^-propagator we have 



ty (fc)^(-fc)> = -^_ . (3.6a) 



In Coulomb gauge, V • A = 0, one finds for the gauge 
field propagator 



(Mk) A 3 (-k)) = 4tt M Sij J ikj , (3.6b) 



where ki denotes the components of the unit vector 
k = k/k. (...) denotes an average with respect to the 
Gaussian part of the action, Eq. (|2.3|) . The vertices as 
given by the action, Eq. (|2.3|) . are shown graphically in 
Fig. [21 and the one-loop diagrams that renormalize the 




FIG. 2: Interaction vertices as given by the action, Eq. (|2,3|) . 
Solid and dashed lines denote the electron fields tj> and tp*, 
respectively, and wavy lines denote the gauge field A. The 
dotted lines serve to separate different components of the elec- 
tron fields. 



C. e-expansion about d — 4 

We first perform a momentum-shell renormalization- 
group (RG) analysis of the action, Eq. (|2.3|) , in d = 4 — e 
dimensions. The propagators can be read off the action, 



various coupling constants in the action are shown in Fig. 
O The calculation is now a straightforward generaliza- 
tion of the one for the s-wave case, Ref. [l|. We define 
the scale dimension of a length L to be [L] = —1, and 
exponents rj and t\a by choosing the scale dimensions of 
the fields tp and A to be |/0(a:)] = (d — 2 + rj)/2 and 
[A(a;)] = (d — 2 + t]a)/2, respectively. We find RG recur- 
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then read 



■ 



- . _Q. 



FIG. 3: One- loop diagrams that renormalize the coupling con- 
stants t, c, u, v, q, and ft in Eq. (|2.3[) . 



sion relations 



di 

di 

du 



o 2 (n + 2)u + 4t> 

2-77 t + 3cY/i+- -4 . 

t + c 

(n + 8)u 2 + 8uv + 8v 2 



dv 
~dl 
dc 

di 

dfx 

~dl 

dq 
~dl 



= (e — 2rf) v — 



(t + c) 
I2uv 



— = —r) c — 3 
1 



2 2 

c q n 
t + c ' 
n (3t + c)c 3 q 2 fi 2 
6 (t + c) 4 



(e - r) A ) q- 



(3.7a) 

3c 2 gV, 
(3.7b) 

(3.7c) 
(3.7d) 
(3.7e) 
(3.7f) 



Here / = In 6 with 6 the length rescaling parameter, we 
have redefined 47T/J, — > /z, and we have absorbed a com- 
mon geometric factor in the coupling constants u, v, and 
u. For w = 0, these flow equations reduce to those of Ref. 
Q], as they should. 

We now look for fixed points (FPs) of the Eqs. (|3.7p . 
Since we are interested in superconductors (as opposed 
to superfluids), we are looking for a FP where the charge 
is nonzero. Equation (|3.7fj) immediately yields 



VA 



(3.8a) 



Anticipating a FP value t* of t that is of 0(e), Eq. (|3.7ej) 
then implies (q 2 ^)* = 6e/n + 0(e 2 ). If we choose rj such 
that c is not renormalized, Eq. (|3.7d[) in turn yields 



r)= -18e/n + 0(e 2 ). 



(3.8b) 



We now look for FP values of u and v that are of O(e). 
Equation (|3.7a[) then indeed yields t* = O(e), so the re- 
maining task is to consider Eqs. (|3.7b|[377cl) . Let us define 
x = u* /ec 2 and y = v* /ec 2 . The remaining FP equations 



(n + 8K + 8xy + 8y 2 - 1 + — x + = 0, 



36 \ 108 



n n 



(3.9a) 



30 



ny z + 12xy - I H j y = 0. 

(3.9b) 

This set of two coupled quadratic equations has four 
solutions. Two of these are given by 



y = o, 



and x a solution of 



., , 36 \ 108 
(n + 8 ^ + — x + — y = 0. 



(3.10a) 



(3.10b) 



Equation (|3. 10b[) is the same condition as in the s-wave 
case, Ref. |l|. It has a real positive solution for n > w 
365.9. In this case, i>q = 0, and Mg = c 2 xe + 0(e 2 ) > 
0. We will refer to these as the s-wave FPs. The only 
other real solutions of Eq. (|3.10b[) occur in the unphysical 
region n < —5.9. 
For y ^ 0, we have 



and y a solution of 

n 2 (n 3 + 8n 2 
-2n(n 3 



1 

12 



96n- 
-38n 2 



36 

n 



ny 



(3.11a) 



1152)y 2 

f 24n - 1728)y 



+ {n 3 + 68n 2 + 1008n + 10368) = 0. (3.11b) 

For positive values of n, this equation has real solutions 
only for n > ss 420.9, which leads to two FPs that 
we refer to as the p-wave FPs. A stability analysis shows 
that of the four FPs found, the only stable one is the p- 
wave FP with the larger value of y. Linearizing about his 
FP yields the critical exponent v for the corresponding 
continuous phase transition. For general n, the expres- 
sion is very complicated. In the limit of large n, one 
finds 



1/1/ = 



(l-l/n + 0(l/n 2 ))e + 0(e 2 ). (3.12) 



We conclude that for large n, and close to d = 4, the 
phase transition in p-wave superconductors is continuous, 
as it is in the s-wave case, and that the p-wave case is in 
a different universality class. The analysis also suggests 
that for physical parameter values, n — 6 and d = 3, the 
transition is unlikely to be continuous. Of course, the 
same caveats as in the s-wave case apply with respect to 
the interpretation of this result. As was shown in Ref. 
1, additional information can be obtained by means of 
an expansion in 1/n in d — 3, and we perform such an 
analysis in the next subsection. 
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D. 1/n-expansion in d = 3 

The technique of the 1/n expansion in a fixed dimen- 
sion was developed by Mai^ for an neutral n-component 
vector field, and it was generalized to the presence of a 
gauge field in Ref. [l|. The basic idea is as follows. Con- 
sider the fully dressed or renormalized counterpart G of 
the Gaussian propagator given in Eq. (j3.6a|) . Its inverse 
can be written in terms of a self energy £, 

G _1 (t, k) =t + ck 2 + E(t, k). (3.13) 

At zero wave number, we have 

G- 1 (t,0) = t + E(t,0) =t R , (3.14a) 

where i R is the renormalized counterpart of t. For t ap- 
proaching its critical value t c , it vanishes according to a 
power law 

*ROC(t-tc)T, (3.14b) 

characterized by the critical exponent 7. This implies 
t c = -£(f = t c ,0) = -S(i R = 0,0), which allows to 
determine t c order by order in some perturbative scheme. 
At criticality, Eq. (|3.13| can thus be rewritten 

G~ 1 (t = t c> k) = ck 2 + £(t = t c ,k) - E(i = t c ,0) 
oc k 2 - 11 = k 2 [l - nlnk + 0(r] 2 )] , 

(3.15) 

with the critical exponent r\. Now consider a pertur- 
bative expansion for X with 1/n as the small param- 
eter. Assuming that r\ is small, it thus can be deter- 
mined perturbatively from the wave number dependence 
of E(t = t c , k) — T,(t = t c , 0). To zeroth order in this ex- 
pansion, there is no contribution, so G _1 (i = t c , k) oc fc 2 , 
and r] = 0. To first order in 1/n, there is a contribution, 
as one would expect from the result of the e expansion, 
Eq. (IX8bl) . 

Similarly, 7 can be obtained perturbatively from the 
behavior of the self energy at k = 0. From Eqs. (|3.13[ 
I3.14j) we have 

ta + E(ta = 0,0)- E(t R , 0) = * - t c cx 4 /7 - (3-16) 

To zeroth order in 1/n one finds, in d = 3, £(£r = 0, 0) — 
1/2 

E(£r,0) oc i R . That is, 7 = 2, which is result for the 
spherical model. If we write corrections to this result in 
terms of an exponent £, 1/7= 1/2 — we have 

E(tR = 0,0)-E(t R ,0) cx 4 /2 ~ C = 4 /2 [l-Clni R +0(C 2 )]. 

(3.17) 



£ can just be obtained as the prefactor of a dependence 

1/2 

of the self energy on tJ lni R . 

For s-wave superconductors, this calculation has been 
performed in Ref. [l|. For the current problem, the calcu- 
lation is straightforward. This calculation yields 

V = ^-+0(l/n 2 ), (3.18a) 
6tt z n 

7 = 2 (1 - - + 0(l/n 2 )) . (3.18b) 
\ tt z n J 

All other static critical exponents can be obtained from 

these two by means of scaling relations. In particular, for 

the correlation length exponent v we have 

352X , , n , , . 

v = 1 - — j - + 0{l/n 2 ). (3.18c) 
in n 

To this order, v is positive for n > 352/37T 2 w 11.9, 
which suggests again that the transition is first order for 
the physical case n = 6 in d = 3. 

IV. SUMMARY 

In summary, we have considered the critical behavior 
of a p-wave superconductor in both an e expansion about 
d = 4, and in a 1/n expansion in d = 3, in analogy to the 
analysis of s-wave superconductors in Ref. [l|. We have 
found that the results are qualitatively the same: both 
methods suggest that, for physical parameter values, the 
superconducting transition is fluctuation-induced first or- 
der in nature. For a model with an n-component order 
parameter, the suppression of fluctuations for sufficiently 
large n leads to a continuous transition that is in a differ- 
ent universality class than the corresponding transition 
in the s-wave case. To first order in e = 4 — d, the critical 
value of n that separates the first-order and second-order 
cases is n C)P = 420.9 for the p-wave case, whereas in the 
s-wave case one has n CjS = 365.9. In a 1/n expansion 
in d = 3, the critical n- value in the p-wave case is equal 
to 11.9, compared to 9.72 in the s-wave case. These val- 
ues have to be compared to the physical values n = 2 
and n = 6 in the s-wave and p-wave cases, respectively. 
The caveats related to critical fixed points that are not 
accessible by either perturbative method that have been 
discussed for the s-wave case^£ apply to the p-wave case 
as well. 



V. ACKNOWLEDGMENTS 

This work was supported by the NSF under grant No. 
DMR-05-29966. 



1 B. I. Halperin, T. C. Lubensky, and S.-K. Ma, Phys. Rev. 2 S. Coleman and E. Weinberg, Phys. Rev. D 7, 1888 (1973). 
Lett. 32, 292 (1974). 



6 



3 G. Nounesis, K. I. Blum, M. J. Young, C. W. Garland, and 
R. J. Birgeneau, Phys. Rev. E. 47, 1910 (1993). 

4 C. Dasgupta and B. I. Halperin, Phys. Rev. Lett. 47, 1556 
(1981). 

5 J. Bartholomew, Phys. Rev. B 28, 5378 (1983). 

6 I. Herbut and Z. Tesanovic, Phys. Rev. Lett. 76, 4588 

(1996) . 

7 I. Herbut and Z. Tesanovic, Phys. Rev. Lett. 78, 980 

(1997) . 

8 K. D. Nelson, Z. Q. Mao, Y. Maeno, and Y. Liu, Science 
306, 1151 (2004). 

9 M. Rice, Science 306, 1142 (2004). 

10 K. Machida and T. Ohmi, Phys. Rev. Lett. 86, 850 (2001). 

11 D. Vollhardt and P. Wolfle, The Superfluid Phases of He- 



lium 3 (Taylor & Francis, 1990). 

A. Knigavko, B. Rosenstein, and Y. F. Chen, Phys. Rev. 
B 60, 550 (1999). 

The A-vertex as given by Sa has a zero eigenvalue, and 
hence the A-propagator does not exist. This problem, 
which reflects gauge invariance, is readily solved by adding 
a gauge fixing term to the action that enforces a particular 
choice of gauge, see, e.g., Ref . ITrj. 

J. H. Chen, T. C. Lubensky, and D. R. Nelson, Phys. Rev. 
B 17, 4274 (1978). 

S.-K. Ma, Phys. Rev. A 7, 2172 (1973). 
L. H. Ryder, Quantum Field Theory (Cambridge Univer- 
sity Press, Cambridge, 1985). 



